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Abstract 

Scalar tetraquark states are studied within the diquark-antidiquark picture in a non-relativistic 
approach. We consider two types of confining potentials, a quadratic and a linear one, to which we 
also add spin-spin, isospin-isospin, and spin-isospin interactions. We calculate the masses of the 
scalar diquarks and of the ground-state open and hidden charmed and bottom scalar tetraquarks. 
Our results indicate that the scalar resonances Dp (2400) and D s (2632) have a sizable tetraquark 
amount in their wave function, while, on the other hand, it turns out that the scalar states 
D* s o(2317) and X(3915) should not be considered as being predominantly diquark-antidiquark 
bound states. We also investigate the masses of light scalar diquarks and tetraquarks, which are 
comparable to the measured masses of the light scalar mesons. 

Key words: Scalar tetraquarks, diquarks, confining potential, hyperfine interaction, non- 
relativistic limit. 


1 Introduction 

One of the most important problems in modern hadron physics is to determine the structure and the 
properties of the newly discovered X , Y, Z states as well other enigmatic mesons, such as _D* 0 (2317), 
Dg(2400), 1)^(2460), etc., see e.g. Refs. [H[2 kSIH and refs, therein. These states cannot be accom¬ 
modated within the simple quark-antiquark picture and are therefore of special interest. 

One possibility is to interpret (some) of these enigmatic mesons as tetraquark states where the con¬ 
stituent objects are a diquark and an antidiquark. Namely, although a diquark cannot be a color singlet, 
the attraction between two quarks can be strong, as various approaches based on one-gluon exchange 
processes [5], instantons 0, lattice calculations [7], and quark-diquark models for the nucleon [5] and 
for baryons in general [5] have shown. Thus, the diquark is an important object for the understanding 
of baryon structure and is also potentially important for the understanding of unconventional mesons, 
most notably tetraquarks. In particular, in this work we are interested in scalar diquarks: these are 
‘good diquarks’ in Jaffe’s terminology uni, with vanishing spin and angular momentum and an an¬ 
tisymmetric flavor wave function of the type [q,q'], where q,q' = u,d,s,c,b (a similar antisymmetric 
combination is realized in color space). 

The masses of heavy tetraquarks as diquark-antidiquark bound states were studied in the presence 
of spin-spin interactions in Ref. Ft, [T2] and later in the comprehensive study of Ref. [T3j. The 
masses of tetraquarks were also calculated in a quark model employing a potential derived from the 
AdS/QCD correspondence [l4j . by using a confining interaction and a meson-exchange potential in a 
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non-relativistic approach m, by implementing the Glozman-Riska (flavor-spin) interaction Hamilto¬ 
nian and SU(3) flavor symmetry breaking [161 . and in the framework of a non-relativistic potential 
model which includes a three-body quark interaction Ell¬ 
in this paper we continue along these lines and calculate masses of (hidden and open) charmed and 
bottom ground-state scalar tetraquarks using two potential models in the non-relativistic limit. As 
a four-body system, a tetraquark state is quite different from a conventional qq meson and we solve 
the problem in a two-step procedure: first, we use a quark-quark interaction Hamiltonian in order to 
obtain the mass of a constituent ‘good diquark’ of the type [■ q,q ']. Second, we regard the diquarks 
as point-like objects and use a diquark-antidiquark interaction Hamiltonian in order to obtain the 
tetraquark masses. In both steps we solve the two-body Schrodinger equation by performing a Taylor 
expansion pS (T9l (20] or by using a variational method. We compare the values of the heavy tetraquark 
masses with the values obtained in previous works and discuss some possible experimental candidates. 
Finally, we focus on the light scalar mesons /o(500), A'q(800), /o(980), and ao(980). These states 
have been, and still are, in the center of a vivid debate concerning their nature: there is now a 
consensus that they are not predominantly quark-antiquark objects [2lj . but that they emerge either 
as dynamically generated molecular-type states, see e.g. [221 [231 EH [251 and refs, therein, and/or 
as tetraquark states as proposed some decades ago by Jaffe (27, [28] and further investigated in Refs. 
[HI US El EH EH ESI E3 [Ml 1301 [33 [33 133113H133113S] • (Note that the quark-antiquark states appear 
in the spectrum but are heavier, since they lie above 1 GeV [37||38]). We apply the very same two- 
step approach described above for a system made of two light diquarks. We evaluate the masses of 
light scalar diquarks and tetraquarks and investigate to what extent the light scalar resonances can be 
described as scalar tetraquark objects. 

The paper is organized as follows. In Sec. 2 we introduce the two potential models and present the 
methods to solve the Schrddinger equation in the presence of hyperfine interactions. Our predictions 
for diquarks and scalar tetraquark masses obtained in the two models are presented and discussed in 
Sec. 3. Finally, a summary and discussion are presented in Sec. 4. 


2 The models 

2.1 The Hamiltonian 

The interaction Hamiltonian for the quark-quark interaction leading to the formation of diquarks is 
given by 

H™(x) = V**(x)+H% lp , ( 1 ) 

where the potential V qq {x) consists of three parts: 

V qq (x) = V conf (x)-^-C . (2) 

The first term V con f(x) is a confining potential (see the next subsections), the second term — t/x is a 
Coulomb-like potential due to one-gluon exchange processes, and C is a constant. The variable x is 
the relative quark-quark coordinate. The quantity H qq p is the hyperfine interaction given by: 

H hyp {x) = H s (x) + H^x) + H S i(x) , (3) 

where Hs(x), Hj(x), and Hsi(x) are spin-spin, isospin-isospin, and spin-isospin interactions, respec¬ 
tively. They read explicitly [Ml 33 EE3 S3 S3i EH]: 

“j>(-/j) ( * i -»2 ) - < 4 > 

h,=a, (^) exp (“^) (fi ' , " i> ' (5> 
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Hsi = A, 


si 


( 7 ^ 7 ) eip (“^ 7 ) (ri ' a)<5 ' f2) ’ 


1 , — 1 L j' 1 o I V^-L ^z/\, t '± ^z ) 3 (6) 

\V 7Ta siJ 

where s, and ti are the spin and isospin operators of the i-t-h quark, respectively, while Ay. and 07,, with 
k = S,I, SI are constants. Note that the operator t z has eigenvalue +5 for the u quark, — \ for the d 
quark, and zero for all other quark flavors. Following Refs. uni mi mm mi, the spatial dependence 
of the hyperfine interaction terms is not modelled by a Dirac 6 function, but by a smooth Gaussian 
function. The hyperfine Hamiltonian is treated as a perturbation which slightly modifies the energy 
levels. 

Next, we turn to the diquark-antidiquark potential. First, we recall that the one-gluon exchange 
potential is such that the quark-antiquark potential and quark-quark potentials are related by V qq = 
2 V qq (this is due to the product of Gell-Mann matrices \ • A j, for details see Refs. [45l ]46) WT \ (48 | [49]). 
When turning to the interaction between a good diquark and a good antidiquark, we assume the 
same form as for a quark-antiquark pair uni- Thus, taking into account the factor 2, we get for a 
diquark-antidiquark system: 

H DD (x)=V DD (x) + H° D p , (7) 

where the potential V DD {x) reads 


V DD (x) = 2V conf (x) - ^ - C . 


2 T 


( 8 ) 


The variable x is now the relative diquark-antidiquark coordinate and Hhyp has the same form as 
Hhyp in Eq. |3]). When applied to (good) diquarks, the isospin operator t z has eigenvalue for the 
diquark [it, g] (with q = s,c,b ), — ^ for the diquark [d,q] (with q = s,c,b ), and zero for the diquarks 
[it, d] and [g, g'] (with g, q' = s, c, b). 


2.2 Quadratic confinement 

In this work, we consider both quadratic and linear potentials in order to model confining interactions. 
First, we study the confining potential in Eq. 0 between two quarks as given by (model 1 ) 

bcon/(-r) — , (9) 


where a is a positive constant. Since the potential is assumed to depend on x only, one can factor out 
the angular part of the two-body wave function. The remaining radial part of the wave function for 
the two-body problem with the unperturbed potential V qq (x) is then determined by the Schrodinger 
equation 


' d 2 2d 
dx 2 x dx 


1(1 + 1 ) 


M x ) 


-2 m[Ei - V qq (x)]4>i(x) , 


( 10 ) 


where 'tpi(x) is the radial wave function, l is the angular quantum number, and m is the reduced mass 
of the two-body system, 


17111712 
717-1 + 1712 


( 11 ) 


with mi and 1712 being the constituent quark (and, subsequently, diquark) masses. Now we solve the 
radial Schrodinger equation for the two-body interaction potential |2]). The transformation 


ipi(x) =x 1 yi (x) 


( 12 ) 


reduces Eq. (flOl) to the form 


d 2 

+ 


e; — 2max 2 + 


2 tot 1(1 + 1) 


<Pi(x) = 0 . 


(13) 
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The radial wave function ipi ( x ) is a solution of the reduced Schrodinger equation for the wave function 
of two identical particles with mass m and interaction potential ©, where 

ei = 2m(E l + C ) . (14) 


The effective potential Ui(x) reads 


, 2 mr 1(1 + 1) 


UAx) = 2 max 2 -b „ 

x x z 

In order to solve Eq. m, we perform a Taylor expansion of Ui (x) around x = xi , 

Ui(x) « Ui(xi) + nf{x - xi) 2 , 
where Xi is such that dUi{x)/dx\ x — xi = 0 and 

2 _ 1 d 2 Ui(x ) 


Of = - 

L o 


2 ! dx 2 


Substituting Eq. © for Ui, Eq. ©, into Eq. m we find 


(15) 


(16) 


(17) 


-j^ipi(x) -flf(x- Xifipi(x) = - [si - Ui(xi)\ ifii(x) 


(18) 


which is the well-known equation for a one-dimensional harmonic oscillator. Namely, for a particle 
with mass m, oscillation frequency u/, energy eigenvalues s' = f n H—) tux ', and spatial wave function 


4>(x), the one-dimensional harmonic oscillator equation reads: 


- n 2 


-<t>{x) = ( x ) ■ 


(19) 


We consider here the ground state of the scalar diquarks (l = n = 0). In this way, upon a comparison 
of Eq. (fl8l) with Eq. (fl9l) . we have: 


mix' 2 me' 

“o = —v— , £o — Uo(x 0 ) = 2 ■ 

Finally, the ground-state energy eigenvalue Eg is obtained using Eq. 

1 


E 0m = -C 


2m 


[C/o(a^o) + n 0 ] 


( 20 ) 


( 21 ) 


with the corresponding ground-state wave function 


ip 0 {x) 



where the constant term in front is a normalization constant. 


( 22 ) 


The very same mathematical problem needs to be solved for the diquark-antidiquark state by treating 
(anti)diquarks as point particles under the influence of the potential ©. The energy eigenvalue E 0 d q 
of the tetraquark ground state n = l = 0 is then calculated in the same way. 
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2.3 Linear confinement 

We also model confinement via a linearly rising potential (model 2): 


Vconf {%) — UX ■ (23) 

The potential (0 is now the well-known Cornell potential. Similarly to the potential of model 1, we 
can factorize the angular part of the Schrodinger equation. Upon substituting the potential (l23l) into 
Eq. (flOl) and using the transformation (fl^l) we obtain: 


dx 2 


<M X ) + 


e/ — 2 max + 


2 rriT l{l + 1) 


ipi(x) = 0 . 


(24) 


We use a variational method to solve the Schrodinger equation for the case l = 0 using the normalized 
test function 

¥>o( x ) = ]j~^ xe ~ p2x * ’ ( 25 ) 

where p is the variational parameter. By minimization of the energy of the system, we calculate the 
energy and also the wave function of the system (for further details of this approach, see Ref. 142]). 
Also in this case, the approach can be easily extended to the calculation of the ground-state energies 
of diquark-antidiquark objects. 


3 Diquark and Tetraquark Masses 

In this section we present the results for diquark and tetraquark masses. We first focus on diquarks 
and then on the corresponding tetraquarks containing at least one heavy quark. Finally, we turn to 
light diquarks and tetraquarks. 

3.1 Diquarks 

The diquark masses obtain a contribution from the constituent quark masses as well as from the 
confining and the spin-isospin-dependent interactions: 

A/ diquark — niqi + rriq 2 + Fo,?? + {Hhyp) , (26) 

where m q i is the mass of i-th quark and Eo, qq is the ground-state energy calculated in the previous 
section. The first-order energy correction from the non-confining potential {Hhyp) is calculated using 
the unperturbed wave function obtained in Secs. 2.2 and 2.3. 

{Hhyp) = d X (p Hhyp ■ (27) 

For the numerical evaluation, we use for model 1 the light and heavy quark masses and the parameter 
r from Ref. 145] , while the parameter a and the hyperfine potential parameters are taken from Refs. 
mm- In model 2, the light and heavy quark masses are still taken from Ref. > 45| but the potential 
parameters are from Ref. (SDj. In both models the parameter C is obtained by fitting it to the 
experimental mass of the p meson. The parameters of both models are summarized in Tab. 1. 

The scalar diquark masses obtained by models 1 and 2 are shown in Tab. 2 and are compared with 
the theoretical works [HIHJGn]- We note that the predictions of the two models are similar to each 
other as well as to previous theoretical calculations. 
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Table 1: Parameters used in our model. 


Parameter 

Model 1 

Model 2 

os 

2.87 fm 

2.87 fm 

As 

67.4 fm 2 

67.4 fm 2 

OS I 

2.31 fm 

2.31 fm 

As/ 

-106.2 fm 2 

-106.2 fm 2 

Ol 

3.45 fm 

3.45 fm 

A/ 

51.7 fm 2 

51.7 fm 2 

m u = rrid 

277 MeV 

280 MeV 

m s 

553 MeV 

569 MeV 

m c 

1816 MeV 

1840 MeV 

m b 

5206 MeV 

5213 MeV 

a 

0.73 fm" 3 

1.23 fm” 2 

T 

0.424 

0.287 

c 

-2.684 fm” 1 

-2.08 fm” 1 



Table 2 

Diquark 

masses (in 

MeV). 


Diquark 

Model 1 

Model 2 

Ref. [13] 

Ref. [30] 

Ref. [51] 

M 

406 

527 

710 

395 

441 

ks] 

678 

784 

948 

590 

659 

M 

1918 

2012 

1973 

1933 

1980 

[sc] 

2147 

2213 

2091 

- 

2120 

m 

5296 

5371 

5359 

- 

5140 

[sb] 

5523 

5563 

5462 

- 

5210 
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Table 3: Masses of open charmed and bottom tetraquarks (in MeV). 


Tetraquark 

Model 1 

Model 2 

Ref. m 

Ref. m 

Exp. mm 

[cqqq] 

2398 

2426 

2399 

- 

Eg (2400) 

[cqqs\ 

2618 

2600 

2619 

2371 

D s * 0 (2317) , E* 0 (2632) 

[csgs] 

2855 

2798 

2753 

- 

- 

[bqqq] 

5763 

5748 

5758 

- 

- 

[bqqs] 

5980 

5901 

5997 

- 

- 

[6ssg] 

6217 

6103 

6108 

- 

- 


3.2 Heavy scalar tetraquarks 

Once the diquark masses are calculated, we can evaluate the tetraquark masses by following the same 
steps. The explicit expression reads 


^-tetraquark — ni diquark T rH antidiquark T -Eq ,DD 4“ ^ 


h dd \ 

n hyp j 


(28) 


The results for open charmed and bottom tetraquarks are listed in Tab. 3 and compared with other 
theoretical predictions mm and experimental candidates [52l 153]. The masses of the tetraquarks are 
indeed very similar in all theoretical approaches, with the exception of [csgs], which in our case turns 
out to be heavier than in Ref. m Our results show that the scalar resonance Dq (2400) may contain 
a sizable tetraquark amount in its wave function. On the other hand, the resonance E* 0 (2317) is too 
light to be interpreted as a [cqqs\ tetraquark (see also Ref. [54] for a discussion concerning conventional 
quark-antiquark charmed scalar states). Another interesting but still controversial state is the so-called 
-D* 0 (2632) meson observed by SELEX [53] , the mass of which fits well to our theoretical predictions. 
The results for hidden charmed and bottom tetraquarks are listed in Tab. 4 and compared with the 
theoretical predictions of Refs. mmmm- Also here, the theoretical results are compatible with 
each other. 

Quite interestingly, the by now established scalar resonance Xf(3915) turns out to be too heavy to be 
a cqcq state and too light to be a cscs state. It is then compatible with being a conventional % c o(2P) 
quarkonium state. 

3.3 Light scalar tetraquarks 

Finally, we apply our formalism to the calculation of the masses of light tetraquarks. The aim is 
to understand if the resonances /o(500), A'(800), /o(980), and oo(980) contain a sizable tetraquark 
amount or not (for experiments concerning these states see Refs. [58] [591 .60, 61, and for theoretical 
works concerning the tetraquark hypothesis Refs. [27] [2B1 ESSO][33][62]). In this framework, the scalar 
diquarks behave under flavor (and also color) transformations as antiquarks, 

[it, d] -O- s, [d, s]-H-it, [s,it]-H-d, (29) 

therefore one can construct a nonet of tetraquarks where the lightest state is the [itd][itd] and corre¬ 
sponds to /o(500), the second lightest are the kaonic-like states [sq^ud], [sd][wd] (with q = u,d and 
/ = 1/2) to be identified with A'g(800) and, finally, the four tetraquarks ([s<7][s<?]) with / = 0, 1 which 
correspond to /o(980) and ao(980). 

Using the parameters of Tab. 1 and the diquark masses of Tab. 2 generates tetraquark masses which are 
100-200 MeV heavier than the states /o(500), A'g(800), /o(980), and ao(980). In order to investigate 


7 








Table 4: Masses of double-hidden charmed and bottom scalar tetraquarks and masses of open charmed 
and bottom scalar tetraquarks (in MeV). 


Tetraquark 

Model 1 

Model 2 

Ref. [55] [53 

Ref. [H] 

[cqcq] 

3807 

3662 

3812 

3723 

[cqcs] 

4043 

3862 

3922 

- 

[cscs] 

4268 

4050 

4051 

- 

[bqbq\ 

10521 

10044 

10471 

- 

[bqbs] 

10747 

10228 

10572 

- 

[6s6s] 

10973 

10412 

10662 

- 


Tetraquark Model 1 Model 2 Ref. [56] 


[cqbq\ 

7162 

6908 

7177 

[cqbs] 

7399 

7106 

7282 

[csbq] 

7397 

7096 

7294 

[cs&s] 

7623 

7285 

7398 


Table 5: Masses of light tetraquark states (in MeV). 


Resonance 

flavor content 

I(JP) 

Model 1 

Model 2 

Ref. [62] 

Ref. [63] 

Exp. 62] 

/o(500) 

[ltd] [fid] 

0(0+) 

546 

614 

550 

596 

400-550 

K (800) 

[ltd] [sd] 

1/2(0+) 

765 

804 

767 

730 

653-701 

/o(980) 

[us] [ms] + [ds][ds] 

0(0+) 

962 

962 

984 

992 

970-990 

ao(980) 

[us] [us] — [ds][ds] 

1(0-) 

984 

984 

984 

992 

983.5-985.9 


whether a better agreement is possible, we re-fit the parameter C of Eq. m for both models, using 
the well-known tetraquark state a o (980) as an input and obtain the light and strange diquark masses 
as following: a) model 1 [ C = —3.507 fm -1 ]: M[ qq j = 244 MeV, M [ qs ] = 515 MeV; b) model 2 [ 
C = —3.045 fm -1 ]: M[ qq j = 330 MeV, M[ ?s ] = 592 MeV. Santopinto and Galata 162! have considered 
a diquark-antidiquark picture of the light scalar tetraquarks in the non-relativistic limit, where the 
masses of the scalar diquarks were obtained as M[ qq ] = 275 MeV, M[ sq j = 492 MeV. 

Using the new diquark masses, we obtain the masses of the light scalar tetraquark nonet listed in 
Tab. 5. Our predictions for the masses of the light scalar tetraquarks are in good agreement with the 
experimental data and also with the results obtained in Refs. [62] and [63j . In our model, a small 
difference between the masses of ao(980) and /o(980) arises from the isospin-dependent hyperfine 
interaction. 

In the context of light scalar states it should be stressed that the role of loop corrections to the self¬ 
energy is surely non-negligible for the masses of these states [2H [231HH [25] . Namely, light scalars have 
a strong coupling to pseudoscalar mesons and a diquark-antidiquark configuration can easily transform 
into a meson-meson one. Moreover, our approach is non-relativistic, thus its application to the light 
scalar sector must be treated with care. Yet, our study shows once more that the light scalar mesons 
are not simple quark-antiquark states but may have a sizable four-quark component. In conclusion, 
we mention that light scalar mesons also play an important role at nonzero temperature [35] and at 
nonzero density [64]. 














4 Summary 

In this work we have calculated the masses of the ground-state heavy and light scalar tetraquarks in 
the framework of a non-relativistic approach with two types of confining potentials, a quadratically 
and a linearly rising one, as well as (iso)spin-(iso)spin interactions. The results for the scalar diquarks 
are shown in Tab. 2, while the heavy scalar tetraquarks are summarized in Tabs. 3 and 4. The results 
of both models are compatible with each other, showing only a mild influence of the particular form 
of the confining potential. Moreover, the results are in agreement, apart from a few exceptions, with 
previous theoretical calculations of Refs. mmm- 

Our results for the masses show that the resonance D* 0 ( 2317) is too light to be predominantly a 
tetraquark state of the type cqqs, while the hidden charmed state X(3915) is too heavy to be cqcq and 
too light to be cscs. On the other hand, the state Ug(2400) and the putative D* 0 (2632) can contain 
an important tetraquark component in their flavor wave function (cqqq and cqqs , respectively). In 
addition to already existing experimental candidates, we also made predictions for the masses of scalar 
tetraquark states which can be discovered in the future (see Tabs. 3 and 4). Namely, some of the X , Y, 
and Z states could turn out to be scalar objects. Finally, we have also studied the light scalar sector 
of QCD and found that the masses of light scalar mesons /o(500), Jig (800), /o(980), and ao(980) can 
be described well in the tetraquark picture (see Tab. 5). 

In this work the masses of the tetraquarks are calculated by using a static approach. Mass shifts take 
place as soon as interactions and quantum fluctuations are taken into account. These modifications 
are typically small for hadrons which are (i) narrow and (ii) are far from any decay threshold. For 
what concerns point (i), the ratio T /(M — E t h), where T is the decay width, M the mass of an hadron, 
and E t h the lowest decay threshold of the state, is an important quantity to estimate the role of loops 
[26] . This ratio is indeed large for the light scalar mesons /o(500) and iFg(800) (see the discussion 
Sec. 3.3), thus loop corrections are surely also an important element towards their understanding. 
Even when this ratio is relatively small (as it usually is for mesons in the cliarmonium region), one 
should then consider point (ii): namely, when M is close to one of the decay threshold (not necessarily 
the lowest), distortion of the spectral functions, mass shifts and sizable meson-meson amounts in the 
wave-function of the unstable meson may occur. This is the case of the light scalar mesons /o(980) and 
do(980): both of them are fairly distant from the lowest threshold [tttt and nr] respectively), but very 
close to the KK threshold. Similarly, the state D* 0 (2317) is pretty close to the DK threshold. More 
in general, many of the newly discovered X , Y , and Z resonances are close to one of their intermediate 
threshold, thus care is definitely needed since the role of loops can be very important. 

In view of this discussion, it must be stressed that also the calculation of decay widths should be 
performed in the future. Namely, it is possible that some of the predicted tetraquark states are, due 
to a ‘fall apart’ decay mechanics, too wide to be measured and that therefore will never be seen in 
experiments. This possibility would explain why only some of the (many possible) tetraquark states 
are actually detectable, that is when the energy threshold of the main decay channel is not too far from 
the mass of the tetraquark state, in such a way that the kinematic suppression balances the large decay 
amplitude. Indeed, this pattern takes place for the light scalar mesons, where /o(500) and /\g(800) 
are very broad, while / o (980) and ao(980) are narrow due to the nearby kaon-antikaon threshold. 
Another (indeed related) improvement is to go beyond the two-step calculations performed in this 
work. Surely, it is much easier to solve two two-body problems than a four-body problem, but a 
general feature of our model (as well as of other tetraquark models) is that the diquarks have a 
dimension which is comparable to that of the tetraquark (about 1 frn). In this respect, there are also 
strong quark-antiquark correlations within the tetraquark, because the diquarks cannot be considered 
as point-like objects. Just as mentioned above, the interchange of a diquark-antidiquark ( qq){qq) 
bound state with a more molecular-like quark-antiquark {qq){qq) surely takes place (and is related to 
the decay of the tetraquark in ordinary mesons). Thus, the view of a pure diquark-antidiquark bound 
state serves as a simple (albeit useful) approximation of the problem, but in the future one should also 
go beyond it and solve a (relativistic) four-body problem. 

In addition to the listed needed improvements, we also mention that our approach can be extended 
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to other quantum numbers as well, thus being potentially interesting to investigate further up to now 
not yet understood mesons. 
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